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ABSTRACT 

We present the results from our deep (16x12 hour) Westerbork Synthesis Radio Telescope 
(WSRT) observations of the approximately 7 square degree Bootes Deep Field, centered at 
14 h 32 m 05!75, 34°16'47'/5 (J2000). Our survey consists of 42 discrete pointings, with enough 
overlap to ensure a uniform sensitivity across the entire field, with a limiting sensitivity of 28/iJy 
(larms)- The catalog contains 3172 distinct sources, of which 316 are resolved by the 13" x 27" 
beam. The Bootes field is part of the optical/near- infrared imaging and spectroscopy survey 
effort conducted at various institutions. The combination of these data sets, and the deep na- 
ture of the radio observations will allow unique studies of a large range of topics including the 
redshift evolution of the luminosity function of radio sources, the K-z relation and the clustering 
environment of radio galaxies, the radio / far-infrared correlation for distant starbursts, and the 
nature of obscured radio loud AGN. 

Subject headings: surveys — radio continuum: galaxies — catalogs 



1. Introduction 

One of the main goals of radio astronomy is to 
fully understand the physics of the population of 
extragalactic radio sources (RSs). Issues include 
the onset and demise of the radio activity and re- 
lated starbursts, the influence of the environment 
on the characteristics of the RSs and the appear- 
ance of the first RSs and their relation to the for- 
mation of galaxies, massive black holes and the 
reionization of the universe. 

Detailed investigation of complete samples of 
bright RSs with redshift information have been 
carried out over the last decades (e.g., 3CRR: 
Laing, Riley & Longair 1983; 6CE: Eales 1985, 
Eales et al. 1997, Rawlings, Eales & Lacy 2001) 
and led to many interesting discoveries. For exam- 
ple, it is now well established that the comoving 
number density of z « 1 powerful RSs is about two 
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orders of magnitude larger than it is locally (e.g., 
Longair 1966, Dunlop & Peacock 1990). Another 
example is that the environment of RSs changes 
with redshift, with bright RSs at higher redshifts 
located in denser environments than locally (e.g., 
Best, Longair & Rottgering 1998). 

Through the selection of RSs that are bright 
and have very steep radio spectra (< —1.3), more 
than 150 powerful galaxies with 2 < z < 5.2 
have been found (De Breuck et al. 2000). The 
large starformation rates (Dey et al. 1997) and 
extremely clumpy optical/IR morphologies (Pcn- 
terrici et al. 1999) provide strong evidence that 
these galaxies are massive galaxies close to the 
epoch of formation. Powerful radio emission is 
most likely caused by accretion onto massive black 
holes (M BH > 10 9 M ; McLure et al. 1999, Laor 
2000), indicating that such massive black holes 
formed alongside or possibly before the formation 
epoch of their host galaxies (e.g., Kauffmann & 
Haehnelt 2000). Recent VLT observations have 
revealed the existence of a large scale structure 
of Ly-a emitting galaxies around the radio galaxy 
1138—262 (z = 2.2), reinforcing the idea that high 
redshift radio galaxies (HzRGs) can be used as 
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Fig. 1. — Composition of the radio source pop- 
ulation as function of 325 MHz flux. The frac- 
tion of flat spectrum sources increases dramati- 
cally with decreasing flux, whereas the steep spec- 
trum component (dominated by AGN) drops pre- 
cipitously. Data are from the source overlap be- 
tween the WENSS and NVSS surveys, at 325 and 
f400 MHz respectively. The errorbars represent 
Poissonian errors, and may be smaller than the 
symbol size in some cases. Note that WENSS sam- 
ple incompleteness sets in around 20 mJy, affecting 
the steep source count most. 

tracers of proto-clusters of galaxies (Pentericci et 
al. 2000). 

At very faint flux densities (e.g., a few tens of 
fijy at 1.4 GHz) the radio source counts are dom- 
inated by the z < 1 starburst population (e.g., 
Richards et al. 1999). Long observations with 
the VLA + MERLIN (Richards 2000; Muxlow et 
al. 1999) and WSRT (Garrett et al. 2000) reach 
such faint levels, and have enabled important con- 
straints to be placed upon the redshifts and nature 
of distant starburst galaxies. 

1.1. Survey Rationale 

Our WSRT survey reaches a la detection 
threshold of 28 /iJy at 1.4 GHz, within a fac- 
tor of 2-3 of the deepest radio observations car- 
ried out so far (cf. Windhorst et al. 1999, 
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Fig. 2. — Expected redshift distribution of the ra- 
dio source population for our survey, based on 
the luminosity functions of Dunlop & Peacock 
(1990, PLE model), and Hopkins et al. (1998, the 
Phoenix survey). 



Richards 2000). The surveyed area is, however, 
large enough (about 7 square degrees) to yield 
enough sources to not be severely affected by low 
number statistics for the RS populations under 
scrutiny. As can be seen in Fig. 1, the compo- 
sition of the radio population changes dramati- 
cally towards lower flux density limits. The higher 
flux levels are dominated by RSs with steep spec- 
tra (< —0.5), and a cross-over to flat-spectrum 
sources appears to occur at the 10-100 mJy (at 
325 MHz) level. The former are mostly identified 
with (powerful) radio sources residing in massive 
ellipticals (e.g., Eales et al. 1997), whereas the 
latter can be tied to a population of starform- 
ing late-type galaxies, especially towards sub-mJy 
1.4 GHz flux density levels (e.g., Windhorst 1999, 
Richards et al. 1999). Therefore, deeper observa- 
tions will not only increase the number of detected 
sources, but it will also provide a better handle 
on the relative makeup of the radio population at 
/iJy flux density levels. 

The radio source population properties do not 
only change as a function of flux density, they also 
vary with redshift. Using radio luminosity func- 
tions for AGN (Dunlop & Peacock 1990) and star- 
bursting populations (Hopkins et al. 1998), we can 
model the expected source counts for a given red- 
shift (and limiting flux density). This is plotted 
in Fig. 2 for a limiting flux density of 140/xJy. Be- 
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yond a rcdshift of z ~ 0.3 the number counts for 
our survey are expected to be dominated by the 
AGN population. The combination of information 
presented in Figs. 1 and 2 makes it clear how de- 
pendent on limiting flux density a perceived radio 
source population is. Indeed, for a limiting flux 
density of 1 mjy (at 1.4 GHz) the population is 
dominated by AGN-type sources at all redshifts 
(Hopkins et al. 1998). 

The survey can detect radio sources at the FR I 
/ FR II break level (10 25 W Hz" 1 at 1400 MHz, 
Owen & White 1991) out to a redshift 4 of z ~ 
3. Fainter FR I type sources with a mean radio 
power of 10 24 W Hz -1 drop out of the sample 
around z w 1, and the much fainter starforming 
systems will not be detected beyond z <~ 0.15 at 
the 10 22 W Hz -1 level. The local starbursting 
system M82 has, for comparison, a radio spectral 
power of 10 21 " W Hz" 1 , given its 1400 MHz flux 
density of 8.363 Jy (White & Becker 1992) and 
the adopted cosmology. Thus the survey provides 
ample data for evolutionary studies of radio loud 
systems out to at least a redshift of 1, whereas a 
more complete census of radio sources down to the 
spiral / starburst level has to be limited to sources 
within z ~ 0.1. 

1.2. Relation to non-Radio Surveys 

Key ingredients for follow-up studies are opti- 
cal / near-infrared identifications and redshift in- 
formation for at least a large fraction of the radio 
sources. For this purpose a number of other sur- 
veys are either being carried out or are planned 
for the same part of the sky. These include: 

The NOAO Deep Wide-Field Survey (Pi's Jan- 
nuzi and Dey). This survey consists of a North- 
ern and a Southern part, with the former field lo- 
cated in Bootes (near the North Galactic Pole), 
covering a 3 x 3 degree region, and the latter lo- 
cated in a 2 x 4.5 degree equatorial region in Ce- 
tus. Both fields have been selected for their low 
mid- to far-infrared cirrus emission, their low H I 
column densities and the availability of high reso- 
lution (~ 5") VLA - FIRST survey radio data. Of 
particular interest to our program is the Bootes 
field, which will be imaged to a limiting surface 
brightness of about 28 th magnitude per square 

4 We adopt H D =65 km s" 1 Mpc -1 , Cl M = 0.3, and f2 A = 0.7 
throughout this paper. 



arcsecond in B, V, and i?, and down to ~ 22 
magnitude per square arcsec in J, H, and K. 
These detection limits will permit the optical and 
near-infrared study of faint, sub-L* galaxies out 
to redshifts of about unity (an L* galaxy will have 
K = 19.82 at z = 1, based on an absolute mag- 
nitude of K = -24.44 ± 0.06, cf. Kochanck et al. 
2001). The typical host galaxies of luminous radio 
sources, with masses well in excess of L* galaxies, 
can be detected out to very large redshifts (based 
on the K — z diagram, e.g., Jarvis et al. 2001, Dc 
Breuck et al. 2002). Given our radio source pop- 
ulation mix of powerful radio sources associated 
with intrinsically bright galaxies at high rcdshift 
and less luminous starforming systems at much 
lower redshifts, we expect to be able to detect opti- 
cal / near-infrared counterparts for most of them. 
The NOAO survey limits are well matched to our 
expected counterpart population. 

The IRAC Shallow Survey (PI. Eiscnhardt). 
The Bootes field will be covered by SIRTF's In- 
fraRed Array Camera (IRAC) in four IR bands 
ranging from 3.6yum to 8^m. Coverage towards the 
longer IR bands up to 160yum will be provided by 
the Multiband Imaging Photometer (MIPS), and 
some spectroscopy by the InfraRed Spectrograph 
(IRS, PI in both cases is J. Houck). 

The NOAO and SIRTF wide field surveys are 
aimed to study, among other things, (i) the evo- 
lution of large-scale structure from z ~ 1 — 4, (ii) 
the formation and evolution of ellipticals and star- 
forming galaxies, and (iii) the detection of very 
distant (z > 4) young galaxies and quasars. The 
SIRTF IRAC and MIPS observations will also de- 
tect star-forming galaxies at mid- to far-infrared 
wavelengths. It is this multi-wavelength aspect of 
the project, covering a large fraction of the elec- 
tromagnetic spectrum (2 radio frequencies, several 
optical and near-infrared bands, and the mid- to 
far-infrared space based SIRTF observations) that 
distinguishes this effort from other deep radio- 
optical/near- infrared surveys like the Phoenix sur- 
vey, (Hopkins et al. 1998; Georgakakis et al. 1999), 
and the Australia Telescope ESO Slice Project 
(ATESP) (Prandoni et al. 2000a,b, 2001). 

1.3. Relation to other Radio Surveys 

The Bootes field has been covered by previous 
radio surveys, most notably by the Westerbork 
Northern Sky Survey (WENSS, Rengelink et al. 
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Fig. 3. — Comparison between the various radio 
surveys covering the Bootes field. The triangles 
represent literature surveys, and the stars are our 
Bootes WSRT and VLA surveys. Overplotted are 
representative radio spectra with varying spectral 
indices: a = —0.25, —0.75, and — 1.25, for flat-, 
normal-, and steep-spectrum radio sources. Any 
given radio source in the VLA-BF survey with a 
slopes shallower than —1.25 should also be present 
in WSRT-BF. On the other hand, of the sources in 
WSRT-BF fainter than the VLA-BF limit (-0.5 
mJy), only the objects with slopes steeper than 
— 1.25 should be present in VLA-BF. 

1997) at 325 MHz, by the NRAO VLA Sky Sur- 
vey (NVSS, Condon ct al. 1998), and by the Faint 
Images of The Radio Sky at Twenty-cm survey 
(FIRST, Becker et al. 1995), both at 1.4 GHz. 
A comparison between the literature surveys and 
our Bootes surveys (WSRT at 1.4 GHz in this pa- 
per, and VLA 325 MHz) is shown in Fig. 3, and 
tabulated in Table 1. The varying survey depths 
and frequencies make combining catalogs to obtain 
spectral index information less than straightfor- 
ward. For instance, combining NVSS and 87GB (at 
5 GHz, Gregory & Condon 1991) only makes sense 
if one is interested in radio sources with strongly 
inverted spectra. For our purpose, since the bulk 
of the radio source population has a spectral in- 
dex of around —0.75 (cf. Fig. 10), a combination 
of surveys like the WENSS and NVSS / FIRST is 



best, as can be inferred from the overplotted com- 
mon radio spectra in Fig. 3. 

However, these surveys do not go deep enough 
to effectively probe the transition in radio source 
population occurring around the 1 mJy level (cf. 
Figs. 1, 9, and, e.g., Windhorst et al. 1999). Our 
WSRT observations do go deep enough, but will 
need low frequency data of matching sensitivity. 
We use the VLA at 325 MHz for this purpose, and 
the data from this program will be described in an 
subsequent paper. However, both the NVSS and 
FIRST survey data have been used to calibrate our 
survey flux densities and positions (cf. Sections 4.1 
and 4.2). 

2. Observations 

The observations were carried out by the West- 
erbork Synthesis Radio Telescope (WSRT) operat- 
ing at 1.380 GHz. The WSRT consists of 14 25m 
telescopes, arranged in a 2.7km East- West con- 
figuration. As the back-end, we used the Digital 
Continuum Backend (DCB) with 8 sub-bands of 
10 MHz bandwidth each. The smallest baseline 
(9— A) was set to 54m to limit shadowing at the 
expense of a reduction in large spatial structure 
sensitivity (~ 800" for this minimum baseline and 
frequency). 

2.1. Field Layout and Instrumental Setup 

We designed a survey layout consisting of 42 
discrete pointings. The separation between the 
grid points was chosen to be 60% of the FWHM of 
the primary beam. Given the used tiling and the 
known attenuation of the beam, a more or less uni- 
form noise-background is obtained with this spac- 
ing (85% of the survey area has a local rms within 
\/2<7 median , cf. Sect 4.3 and Fig. 8). This strategy 
is similar to the one used for the Australia Tele- 
scope Compact Array (ATCA) ATESP survey (see 
Prandoni et al. 2000a for a detailed description). 
The total number of pointings (42) was dictated 
by the need to cover the 325 MHz VLA primary 
beam with a uniform sensitivity (noise) level. 

During each observing block of 12 hours, the 
telescopes were continuously cycled between 3 in- 
dividual grid positions. This "mosaicing" mode 
(Kolkman 1993) allows multiple fields to be ob- 
served, while still retaining a 12 hour uu-coverage 
for the fields individually (albeit sampled non- 
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continuously). The basic integration time for the 
observations is 10 seconds. A typical observing 
cycle can be broken down to a 10s slew-time be- 
tween the grid positions, and 5 x 10s on-source 
time. Even though the net slewing time is less 
than 10s, some extra time is needed for the ar- 
ray to settle itself after the move. The observ- 
ing efficiency with this scheme is therefore 83.3% 
resulting in 3 x 200 m net observing time per 12 
hour cycle. Given the total allocated time for this 
project (192 hours), we used the remaining two 
12 hour blocks to cycle through all 42 positions. 
In this setup, each position was revisited every 42 
minutes (instead of every 3 minutes), resulting in 
a rather sparse uw-sampling and less than 15 m on- 
source time per block. 

Each 12 hour block was sandwiched between 
two phase and polarization calibrators, typically 
3C 286 and 3C 147, more than adequate given the 
system stability. A log of the observations can be 
obtained from the ftp-site (cf. Sect 4). 

3. Reduction 

The mosaic was reduced, calibrated, and as- 
sembled using the MIRIAD (Sault et al. 1995) 
software package. The data was typically of high 
quality, and only a few percent had to be flagged. 
Usually the bad data was concentrated in channel 
5 (of the 8), which, around 1380 MHz, is the fre- 
quency band most affected by interference due to 
the Global Positioning System. Every field point- 
ing was mapped using a multi frequency synthe- 
sis approach, where the measurements of the 8 
bands individually are gridded simultaneously in 
the uv plane. This significantly reduces bandwidth 
smearing problems. Then a three step iterative 
phase self-calibration cycle was used, using typ- 
ically around 100000 clean iterations. The clean 
was done down to the 3 sigma level, so as not to in- 
corporate too much flux in faint sources that does 
not belong there. In the few fields with strong 
sources present we performed amplitude self cali- 
bration as well, in all other cases this did not im- 
prove the final outcome. In the fourth, and final 
cycle, spectral index effects on beam-shape were 
taken into account (Sault & Conway 1999). The 
final maps improved significantly by correcting for 
the small primary beam shape variations across 
the 8 10 MHz wide frequency channels. 



After all the field pointings were reduced in 
this manner, they were assembled into a final mo- 
saic. This step basically adds up the maps after 
performing the proper primary beam correction. 
Since the dirty beam changes slightly across the 
field, we restored all the data with a fixed syn- 
thesized beam of 13.0" x 27.0", at a position an- 
gle of zero degrees. The mosaic was then further 
mapped onto 4" x 4" pixels, for a total of ~ 2750 2 
pixels. 

4. Results 

We used automated routines for the source ex- 
traction and catalog creation. These were slightly 
modified from their WENSS counterparts, but the 
applied methods are exactly the same, all of which 
are described in detail in Rengelink et al. (1997). 

The software works on rectangular patches of 
sky only, so we tiled the circular overall shape of 
the survey into three rectangular areas, as outlined 
in Fig. 8. All of the low noise areas have been in- 
cluded this way, and only a few parts of the noisy 
edges have not been cataloged. The total cata- 
loged survey area covers 6.68 square degrees. 

Part of the catalog, to illustrate its format, has 
been listed in Table 2. The full version (with 
3172 sources) can be obtained through anonymous 
ftp to ftp://ftp.nfra.nl/pub/Bootes. The com- 
plete mosaic, individual pointing maps, and tables 
with various additional data are available from the 
same address, all of which are described in the 
README file. 

4.1. Flux Accuracy and Error Estimates 

We have compared the flux densities of unre- 
solved sources present in both our uncalibrated 
Bootes and the NVSS catalog. Since the reso- 
lution of NVSS (at 45") is slightly worse than 
ours (13" x 27"), a source which is unresolved in 
our catalog is consequently unresolved in NVSS. 
While we could compare our fluxes to the deeper 
FIRST data, the latter's much higher angular res- 
olution typically resolves point sources in our cat- 
alog, making a direct comparison difficult. The re- 
sults of the comparison are plotted in Fig. 4. It is 
clear from the plot that our uncalibrated fluxes are 
a little too high in comparison to the NVSS fluxes, 
at least for S > 10 mJy. Below these fluxes, the 
NVSS values are systematically too high, and are 
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Fig. 4. — Flux density comparison between 
sources in common to NVSS and WSRT. Only 
unresolved WSRT sources (and hence unresolved 
in NVSS) have been included. The increase in 
NVSS flux density relative to the WSRT close to 
its detection limit is possibly due a combination of 
Malmquist and NVSS clean biases. Errorbars are 
la errors, and may be smaller than the symbol 
size in some cases. 



presumably due to a combination of Malmquist 
and clean biases. This overestimate of NVSS fluxes 
close to their detection limit (<~ 2.5 mJy) is also 
evident in the Condon et al. (1998) comparison of 
NVSS fluxes to deep WSRT (Katgert-Merkelijn et 
al. 1985) flux densities, cf. Fig. 31 in Condon et al. 
(1998). Also, Prandoni et al. (2000b) noticed the 
same effect in comparing their ATCA radio survey 
fluxes to the NVSS values. 

Using flux density weighting, we calculated the 
offset to be 4.5 ± 3.0% too high. We reduced our 
fluxes accordingly (cf. Fig. 5). Following Ren- 
gclink et al. (1997), the relative flux density errors 
can be written as: 



S 



2 / Crms \ ' 



V S J 



(1) 



This equation reflects the two components of the 
measurement error, with C\ due to a constant 
systematic error, and C2 being dependent on the 



Fig. 5. — Relative flux density errors, as function 
of Signal-to-Noise ratio. WSRT fluxes are com- 
pared to NVSS fluxes for unresolved radio sources 
present in both catalogs. The overplotted solid 
lines are expected upper and lower envelopes to 
the flux density fraction, assuming values for C\ 
and C2 (cf. Eqn. 3) of 0.04 and 1.3 respectively. 
The dashed lines are for C\ = 0.0 and C2 = 2.0, 
and illustrate that C\ is most dominant at high 
SNR while C 2 is it at low SNR. 

Signal-to-Noise Ratio (SNR). Ideally, one would 
like to compare our measured (and corrected) 
fluxes to their true values in order to determine the 
constants C\ and C2 which best fit the observed 
>S'wsRT/>S'truc ratio. Unfortunately, we do not have 
such a control data set, instead we use the NVSS 
measurements. If we assume that the NVSS mea- 
surements have a similar error dependence, we can 
define the flux density ratio as: 



1 ± (cl + cf ( £ gf t ) 



2\ 2 



(2) 



with CTwsrt and <r NVSS being the median noise in 
the sky, and 5 tru6 the true value of the source 
flux. The ct's are quoted as 0.45 mJy for the 
NVSS (Condon et al. 1998), and 0.028 mJy for 
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our survey. If we further assume that S WSRT and 
S'nvss are approximately equal to S ttue , and that 
SNRx = Sxjox with X being either WSRT or 
NVSS, we can rewrite Eqn. 2 as: 

w_ i*(c? + Glk) 2 )* (3) 

based on ct NVS s = 16.1 x cr WSRT , which implies 
SNR W 

'srt — 16.1 x SNR NVSS . The results using 
Eqn 3 has been overplotted on Fig. 5 such that 
the maximum (upper) envelope is given by setting 
the ± to 1 + ()/l — (), and the minimum envelope 
by 1 — 0/1 + in the equation. The C parameters 
have been set to 0.04 and 1.3 respectively, identi- 
cal to the values in Rcngelink ct al. (1997). The 
model is most sensitive to the C\ value, which ba- 
sically sets the envelope separation at high SNR. 
The C2 value, which scales the SNR dependence 
is far less constrained. Values of C 2 = 2 (e.g., 
Kaper et al. 1966) are not excluded. Given the 
assumptions and the assumed uncertainties about 
the NVSS errors, we adopt the WENSS values of 
0.04 and 1.3. The quoted flux density errors in the 
final catalog are calculated with these particular 
values. 

4.2. Positional Accuracy 

Optical identifications can only be securely 
made if the radio positions are known accurately. 
The optical source density becomes high enough 
towards fainter magnitudes to effectively have one 
potential counterpart per beam. For instance, the 
mean source separation in the Deeprange /-band 
field survey is about 17" at the 23 rd magnitude 
level (Postman et al. 1998). This separation is ac- 
tually smaller than the WSRT beam size at 1400 
MHz. Good positional matches are therefore es- 
sential. 

We compared the cataloged positions for point 
sources against their FIRST positions and against 
Automatic Plate Measuring (APM) machine iden- 
tifications. The APM facility (in Cambridge, UK) 
catalogs identifications and positions based on 
scanned UK and POSS II Schmidt plates, covering 
currently more than 15 000 square degrees of sky. 

The relative offsets for the individual sources 
are plotted in Fig. 6. It is clear that the FIRST and 
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Fig. 6.— WSRT positional check against FIRST 
and optical APM positions. Only WSRT point 
sources are used. FIRST positions are indicated 
by the solid circles, and APM identifications are 
potted as open circles. The larger scatter in the 
APM correlation is due to the association of ra- 
dio positions with unrelated nearby optical ob- 
jects. The mean offsets (indicated by the crosses) 
are ARA= +0.60", ADEC= -0.17" for FIRST 
and ARA= +0.51", ADEC= +0.02" for the APM 
match. 

APM positions agree rather well with each other 
(indicated by the crosses), but that our positions 
are systematically off in RA. Without more fre- 
quent observations of additional calibrators, which 
would adversely affect our uv coverage, astromet- 
ric accuracy of the WSRT is known not to be better 
than about 0.5" (e.g., Oort & Windhorst 1985), 
consistent with our offset value. We corrected all 
the positions in RA with —0.56", i.e., the mean 
of the FIRST and APM RA offsets. This correc- 
tion corresponds to about 4% of the beam width, 
small but significant enough when accurate posi- 
tional coincidences are needed. 

Analogously to Eqn. 1 for the flux density er- 
rors, we can describe the flux density dependence 
on positional accuracy in the form: 



7 



u i i °j i r~ i i - i i s i i i i i i 

12 3 4 

LOG(SNR) 

Fig. 7.— Positional offsets from FIRST / APM 
positions as a function of the signal-to-noise ra- 
tio (SNR). Overplotted as solid squares are the 
67 percentile values for a given SNR. The bottom 
curve represents the best fitting la error envelope 
(cf. Eqn. 4, with d = 0.44 and C 2 = 5.5. The 
top curve is identical to the one modeled for the 
WENSS survey (Rengelink et al. 1997). 



= (cf + Cf (SNR)~ 2 ) ' (4) 

The absolute distances from the FIRST posi- 
tions have been plotted in Fig. 7 as a function of 
signal-to-noise. Our survey positions have been 
corrected for the RA offset first. Since the FIRST 
resolution is higher than our survey and may lead 
to 2 (or more) FIRST catalog positions for any 
of our positions, we only considered FIRST point 
sources within our survey field. The inclusion of 
resolved sources (cither in FIRST or our survey) 
unnecessarily complicates the comparison. 

In Fig. 7 a clear decrease in positional offset 
with increasing SNR can be seen. To characterize 
this trend, we fitted Eqn. 4 to the 67th percentile 
points (solid squares) in order to get a la posi- 
tional error estimate. The best fitting values for 
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Fig. 8. — Map of the rms noise inside the source 
extraction area (outlined by the boxes) . The noise 
levels and relative surface area are, from inside 
out: 20-30^Jy (60.9%), 30-40^Jy (23.9%), 40- 
80^Jy (13.3%), and 80 - 160/iJy (2.0%). 

the constants are: C\ = 0.44 and C2 = 5.5, both in 
arcseconds. An outer envelope to the offset distri- 
bution is given by C\ = 0.44 and C 2 = 15.0. The 
value for C2 is actually the mean beam size (taken 
to be 20") divided by 1.3, a value identical to the 
one quoted for the WENSS survey (Rengelink et 
al. 1997). We adopt the first set of constants (the 
la equivalents) for our source catalog. 

4.3. Completeness and Reliability 

The background noise in our survey is not uni- 
formly flat, but has a marked upturn towards the 
edges. The tiling was set up in such a way that in 
the interior regions the noise should be flat. This 
can be verified in Fig. 8, which plots the actual 
background noise. The median noise level of the 
inner parts is 28/xJy. The large 30 to 40^iJy "in- 
trusion" at 14 h 36 m , 34 d 00 m is most likely due to 
the somewhat higher noise levels in those four par- 
ticular pointings. 

The survey completeness can be gleaned from 
Fig. 9, in which the differential source counts are 
plotted against flux. The number counts have 
been normalized by the expected number in a Eu- 
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Fig. 9. — Differential source counts as function of flux density S (vertical bars). The counts have been 
normalized to the expected number for a Euclidean universe. The lengths of the vertical bars are set by 
the >/N error-bars. Overplotted are the differential source counts for the NVSS (stars) and FIRST (open 
squares) surveys. Note the onset of incompleteness in NVSS at 12 mJy, in FIRST around 2 mJy, and in our 
survey at about 0.2 mJy. 



clidean universe, given by AT-g c i = K V SZ , with 
the unitless constant K v = 225 (consistent with 
e.g., Oort & Windhorst 1985 and Oort 1987, but 
see Wall (1994) who lists K v = 200, however). 
Our source counts are compared to the ones based 
on the NVSS and FIRST catalogs, which due to 
their much larger survey area extend further to- 
wards higher flux densities. There is good agree- 
ment (within the la errorbars) over the range 10 
- 100 mJy between the surveys. The small devi- 
ation in our differential counts around the 4 mJy 
bin appears to be real and might indicate the pres- 
ence of an overdense region within our field (e.g., 
a cluster). Since our survey field is relatively small 
at 6.68 square degrees, any local overdensity could 
skew the number counts significantly. The NVSS 
and FIRST number counts are not affected by this, 
and serve as a useful baseline. 

All three plotted surveys have different com- 
pleteness limits, and if one marks the first system- 
atic deviation from a low-order polynomial fit as 
the completeness limit, we measure 12 mJy for 
NVSS, 2 mJy for FIRST, and 0.2 mJy for our 
survey. However, since the noise in our survey 
is not constant across the source extraction area, 



but varies within the 30% level over ~90% of the 
survey, the 0.2 mJy value is not strictly correct. 
It represents a mean value over the survey area, 
with completeness levels slightly lower and higher 
in the inner and more outer parts of the survey 
respectively. 

4.4. Source Confusion 

With the 13" x 27" beam size and the faint flux 
density levels reached in this survey, considerable 
source confusion might be present. We therefore 
modeled this by randomly distributing the cata- 
loged source population over the survey area a 
large number of times (10 4 ~ 5 ). Each time the 
likelihood of having close pairings of sources was 
recorded. Given a large enough sampling, a more 
or less accurate estimate of the frequency of oc- 
currence is possible. The results are given in Ta- 
bic 3. Since we used the actual catalog, a strong 
flux dependency is to be expected. In other words, 
having two bright (> 10 mJy) objects very close 
together almost always means they are physically 
associated, whereas two faint (< 0.5 mJy) objects 
with a similar separation are most likely uncon- 
nected. Based on the numbers from Table 3 we 
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can state that sources with a double morphology, 
with flux densities > 3 mJy, and angular separa- 
tions of < 1 arcminute have a 96% chance of being 
true physical doubles. All of the listed resolved 
sources with flux densities > 10 mJy in Table 2 
should therefore be considered single physical en- 
tities and hence accurately classified. 

Assuming the approximately 2-fold increase in 
unrelated object count continues toward fainter 
flux levels, a similar WSRT survey would become 
confusion dominated (i.e. with on average 2 faint 
objects in the synthesized beam) around the 15 
/iJy mark. For this survey, with a 5a limit of 140 
^Jy, on average 10% of the beams are confused. 

4.5. Source Catalog 

The final catalog, which lists every source with 
a flux density over 5 times the local a (correspond- 
ing to 140^iJy in the center of the survey field), 
contains 3172 sources. Roughly 10% of these are 
resolved (316) by the 13" x 27" beam. A com- 
plete break down of source morphology is given 
in Table 4. More complete catalogs with varying 
threshold cr's are available from the ftp site. The 
total number of included sources decreases with in- 
creasing SNR limits: 3172, 2767, 2367, 2061, 1854, 
and 1692 sources for 5, 6, 10er thresholds re- 
spectively. Also, detailed radio maps of the com- 
plete survey area are provided, with the cataloged 
sources clearly indicated. This will allow for a di- 
rect visual assessment whether a particular source 
is to be considered real or not. 

The 73 brightest resolved objects have been 
listed in Table 2, all of which have flux densities 
in excess of 10 mJy. Contour plots for these par- 
ticular sources are presented here in Fig.ll. 

Our survey catalog contains 143 sources which 
are also detected in the WENSS survey. Fig. 10 
plots the 325— 1400MHz spectral index distribu- 
tion of these sources. There appears to be a trend 
for the more luminous sources (> 100 mJy) to have 
slightly steeper spectral indices than the fainter 
part of the sample (< 100 mJy). The actual mean 
values are -0.60 ± 0.31 and -0.81 ± 0.13 for the 
flux bins 10 - 100 and 100 - 1000 mJy respec- 
tively. The quoted errors are the la standard de- 
viations. This overall flattening of the spectral 
index with decreasing flux density levels is con- 
sistent with the data presented in Fig. 1, which 




10 100 
1400 MHz flux [mJy] 



Fig. 10. — Spectral index distribution func- 
tion of 1400MHz flux. The 325MHz data is from 
the WENSS survey. The slanted line represents 
the WENSS sensitivity limited spectral index, us- 
ing a 3a detection threshold of 11 mJy at 325MHz. 
The two solid horizontal lines are the spectral in- 
dex means for the flux density ranges 10 to 100, 
and 100 to 1000 mJy, respectively. The dashed 
boxes outline the la area. 

shows the change in radio source population as a 
function of flux density based on the WENSS and 
NVSS surveys. Unfortunately, we cannot use our 
much deeper survey (compared to NVSS) to ex- 
tend this towards even lower flux densities. To 
the left of the slanted line in Fig. 10 the WENSS 
survey was not deep enough to detect the radio 
sources at 325MHz. We will use our deep VLA 
observations of the Bootes field for this purpose. 

Ten of the radio sources in our catalog have red- 
shifts given in the literature, and these are listed 
in Table 5. Alternative (radio catalog) names for 
some of our objects are given in the footnotes to 
Table 2. 

5. Summary 

We presented the results from our deep WSRT 
observations of the Bootes Deep Field. The sur- 
vey reached a la limiting flux density of 28^iJy 
in the central region, and 3172 sources were de- 
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tected above the 5c level in a 6.68 square degree 
area. The survey is deep enough to sample the 
change in radio source population properties at 
the few mjy level. In combination with our lower 
frequency VLA data, these datasets will provide 
key information pertaining, among other things, 
to the nature and evolution of radio sources, both 
in the local and the high redshift universe. 
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Fig. 11. — Plots of resolved sources with S > 10 mJy. The contour levels are given by (1 + 2 n )a, with 
n = 1, 2, ... and a as the local background noise. Negative contours have the same spacing, and are plotted 
as dotted lines. The cross represents the nominal source center, and the object name and 1400 MHz flux 
densities are given in the upper left and right corners. 
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Table 1 

Comparison between various radio surveys covering the Bootes field 



Survey 


Frequency 


Resolution 


Flux density limit a 


Detections in Bootes Field b 


WENSS 


325 


54" x 99" 


15 


180 


VLA-Bootes 


325 


6" 


0.5 


1200 c 


NVSS 


1400 


45" 


2.5 


438 


FIRST 


1400 


5" 


1.0 


749 


WSRT-Bootes 


1400 


13" x 27" 


0.140 


3172 


87GB 


4850 


222" x 198" 


18 


22 



a 5cr detection limit. Units are in milli-Jansky's. 

b Thc number of radio sources/components within a circular aperture with a radius of 5400", centered on 
14 h 32 m 05!75, 34°16'47"5. Note that this number is both depending on flux limit and resolution. 

c Expected. 



19 



Table 2 

Resolved sources with flux densities >10 mJy 



Source 


RA (J2000) 


Dec (J2000) 


6POS 
["] 


-pa 




F int 
[mJy] 






rms b 
[mJy] 


©ma 
["] 


["] 


PA 

[°] 


LAS d 
["] 


J 14z041 + o40oz0 


14 


25 


41. 


.40 


+34 


58 


26. 


.6 


o 


.4 


A/T 


OfiO 
zOz 
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± 


10.5 


0. 


.102 




QO 

oy 


07 
Z / 


OA Q 
Z4o 


J 14zOU 1 +o4U4z4 


14 


2G 


07. 


,73 


+34 


04 


24. 


.0 





.4 


A/T 
1VI 


97 


73 


± 


1 


.51 


0. 


.048 


111 
111 


A O 

4y 


1 RA 
104 


1 nc 
lUo 


J 14z01U+oooyoO 


14 


26 


10. 


.88 


+33 


39 


36, 


6 


o 


.4 


A/T 
1V1 


zy 


.78 


± 


1 


.19 


0. 


050 


1 KCi 

ioy 


fi7 


1 
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1 K1 
10 / 


J 14z0zU+o44UlZ 


14 


26 


20. 


94 


+34 


40 


12, 


6 


o 


.4 


A/T 


Q7 
O I . 


,01 


± 


3 


.48 


0. 


,034 


OQ1 

zol 


77 


1 QO 
loU 


000 

zzy 


T1 /10£QO_LQ/1/1Q1 Q 
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14 


26 


39. 
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+34 


43 


18, 


.6 


o 


.4 


A/T 
Ivl 


1 £ 
10. 
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± 





.64 


0. 


,031 


1 KC\ 
lOU 


01 

y i 


OO 
ZZ 


1 A Q 
14o 


Ti/io^^n i qqoO/11 
J 14z00U+oozy41 
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26 


■50. 
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41, 


,6 


o 
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1VI 


1 Q 
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± 
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.31 
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1 11 
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1 70 
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1 IK 
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Table 2 — Continued 
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a Flag: S=point source, M=resolved, E=complex. 
b Local sky RMS (in mJy). 

c Apparent angular extent of 2.5cr contour. Sources with ma j/6min ~ 2 and PAw 0° are considered barely resolved. 

d Largest Angular Size. Resolved sources are deconvolved with the beam size, point source sizes are approximated by: 
LAS = e bo am x {0.04 2 + (6.0/SNR) 2 }!, cf. Rcngelink ct al. (1997). 

Alternative names: J142659+341159 = 7C 1412+344; J142739+330750 = 7C 1425+333; J143025+351914 = 
NGC 5656; J143048+333319 = 7C 1428+337; J143052+331320 = 7C 1428+334; J143134+351506 = 7C 1429+354; 
J143317+345108 = 7C 1431+350; J143342+341138 = 7C 1431+344; J143749+345452 = 7C 1435+351 
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Table 3 
Source Confusion Limits. 



Lower flux limit Expected Number of Objects within radius R 
[mJy] R = 30" R = 60" R = 120'' 



0.05 


0.11 


0.43 


1.70 


0.10 


0.10 


0.43 


1.69 


0.20 


0.074 


0.29 


1.18 


0.40 


0.041 


0.17 


0.65 


0.80 


0.025 


0.097 


0.38 


1.60 


0.012 


0.057 


0.23 


3.20 


0.0083 


0.036 


0.14 


6.40 


0.0052 


0.019 


0.073 


12.80 


0.0026 


0.0090 


0.039 


25.60 


0.0012 


0.0055 


0.021 


51.20 


0.0007 


0.0024 


0.0095 


102.40 


0.0002 


0.0005 


0.0031 


204.80 


0.00002 


0.0001 


0.0007 



Note. — The listed counts are the number of unrelated objects within 
a search radius R around a preselected target. For the total source count 
within a radius R, 1 should be added to this count therefore. The lower 
2 flux density bins are affected by the incompleteness of the catalog at 
those levels; otherwise a factor of <~ 2 decrease in expected counts with 
increasing flux density threshold seems to be present. Also note the sur- 
face area factor of 4 in count levels between the columns. Objects within 
the 30" radius are too close to be resolved by the WSRT beam, and would 
mistakenly be classified as a single source. To keep the number of sig- 
nificant digits approximately constant we had to increase the number of 
simulations with increasing flux density threshold. 



Table 4 

Catalog morphology break-down. 



Morphology 


Number 


Density / sq. degree 


Example object 


Unresolved 


2856 


427.5 




Barely resolved 


43 


6.4 


J142851+353029 


Double 


136 


20.4 


J143703+343442 


Triple 


13 


1.9 


J143309+333609 


Asymmetric 


112 


16.7 


J143604+334539 


Complex / Other 


12 


1.8 


J143429+342812 


Total 


3172 


474.9 
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Table 5 

Sources with known redshifts. 



Name 


Alt. Name 


RA (J2000) a 


DEC a 




ID 


z 


V 


f 1400 


J142744+333829 




14 27 44.49 


+33 


38 29.2 


G 




1.237 




14.18 


J142823+331514 


UGC 09284 


14 28 23.42 


+33 


15 14.2 


G: 


Sa? 


0.01386 


14.86 


0.93 


J142932+333038 


CG 0447 


14 29 32.66 


+33 


30 38.4 


G 




0.02638 


17 


0.12 


J142934+352742 


NGC 5646 


14 29 34.07 


+35 


27 42.2 


G: 


SBb 


0.02861 


14.99 


1.37 


J143025+351914 


NGC 5656 


14 30 25.43 


+35 


19 14.6 


G: 


SAab 


0.01051 


12.73 


15.33 


J143119+343803 


CG 0457 


14 31 19.91 


+34 


38 03.9 


G 




0.01440 


17.20 


1.25 


J143125+331349 


CG 0459+0460 


14 31 25.36 


+33 


13 49.9 


G: 


S 


0.02247 


14.6 


4.73 


J143156+333830 


VV 775 


14 31 56.15 


+33 


38 30.1 


G: 


Irr 


0.03373 


16 


5.88 


J143232+340626 


LEDA 099838 


14 32 32.42 


+34 


06 26.3 


G 




0.04264 


17.96 


0.21 


J143518+350709 


CG 0479 


14 35 18.28 


+35 


07 09.2 


G 




0.02847 


14.5 


26.88 



a Radio position. 

b Radio flux in mJy, with a = 0.028 mJy. 
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